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Abstract. Microscopic theory of counting statistics of electrical noise is 
reviewed. We discuss a model of passive charge detector based on current 
fluctuations coupled to a spin, and its relation with the theory of photon 
counting in quantum optics. The statistics of tunneling current and, in 
particular, the properties of the third moment are studied in detail. The 
third moment is shown to be temperature-independent for tunneling in a 
generic many-body system. Then the statistics of mesoscopic transport is 
discussed. We consider applications of the functional determinant formula 
for the generating function of counting distribution to the DC and photo- 
assisted transport, and to mesoscopic pumping. A universal dependence 
of the noise in a mesoscopic pump on the pumping fields is obtained and 
shown to provide a method of measuring the quasiparticle charge in an 
open system without any fitting parameters. 
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1. Introduction 

The measurements performed by optical detectors, such as photon counters, 
are extended in the time domain, which makes them sensitive to tempo- 
ral correlations of photons [1]. It has been known long ago in the theory of 
photodetection [2] that understanding photon counting is essentially a prob- 
lem of many-particle statistics. Similar considerations apply to the electri- 
cal noise measurement, although it differs from photodetection in that the 
electrons, unlike photons, are not destroyed in the process of counting. The 
noise measurement, very much like photodetection, is a sensitive probe of 
temporal correlations between electrons. 
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Fermi correlations in the electron noise were originally studied by Leso- 
vik[3] (see also Ref. [4]) in a point contact, and then by Biittiker [5] in 
multiterminal systems, and by Beenakker and Biittiker [6] in mesoscopic 
conductors. Kane and Fisher proposed to employ the shot noise for detect- 
ing fractional quasiparticles in a Quantum Hall edge system [7] . Subsequent 
theoretical developments are summarized in a recent review [8] . 

Experimental studies of the shot noise, after first measurements in a 
point contact by Reznikov et al. [9] and Kumar et al. [10], focused on the 
quantum Hall regime. The fractional charges e/3 and e/5 were observed [11, 
12, 13] at incompressible Landau level filling (see also recent work on noise 
at intermediate filling [14]). The shot noise in a mesoscopic conductor was 
observed by Steinbach et al. [15] and Schoelkopf et al. [16], who also studied 
noise in photo-assisted phase-coherent mesoscopic transport [17]. 

In this article we discuss counting statistics of electric noise and con- 
sider the probability distribution of charge transmitted in a fixed time inter- 
val [18, 19]. This distribution provides detailed information about current 
fluctuations. The counting statistics have been studied for the DC transport 
of free fermions [18, 20], the photo-assisted transport [21], the parametri- 
cally driven transport [19, 22], and in the mesoscopic regime [23] (also, see a 
review [24] ) . Nazarov developed Keldysh formalism for the counting statis- 
tics problem and applied it to mesoscopic transport in a weak localization 
regime [25] and, together with Bagrets, in a multiterminal geometry [26] . 
Charge doubling due to Andreev scattering in NS junctions was considered 
by Muzykantskii and Khmelnitskii [27], and in mesoscopic NS systems by 
Belzig and Nazarov [28] . Andreev and Kamenev [29] studied the problem of 
mesoscopic pumping in view of the results of Ref. [19]. Taddei and Fazio 
discussed counting statistics of entangled electron sources [30] . Statistics 
of transport in a Coulomb blockade regime was studied by Bagrets and 
Nazarov [31]. Photon statistics was considered by Beenakker and Schome- 
rus [32] and Kindermann et al. [33] . The problem of back influence of a 
charge detector on current fluctuations in the context of counting statistics 
measurement was studied by Nazarov and Kindermann [37] . 

The possibility of measuring counting statistics using a fast charge inte- 
grator scheme was considered recently [34] . From the measured distribution 
all moments of charge fluctuations can be calculated and, conversely, the 
knowledge of all moments is in principle sufficient for recovering the full 
distribution. However, due to the central limit theorem, high moments are 
probably difficult to access experimentally. Therefore recent literature fo- 
cused primarily on the third moment. It was found that the third moment 
obeys a generalized Schottky relation which holds in the tunneling regime 
at both high and low temperature, but involves a temperature-dependent 
Fano factor in the mesoscopic regime [22, 35]. Gutman and Gefen [35] stud- 
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ied the third moment using a sigma model approach, while Nagacv [36] 
demonstrated that all moments are correctly reproduced by an extension 
of the Boltzmann-Langevin kinetic equation. 

In this article, after introducing the counting distribution (Sec. 2), we re- 
view its microscopic definition based on a passive charge detector (Sec.2.1). 
In Sec. 2. 2 we study the statistics of tunneling in a generic many-body sys- 
tem. From the microscopic approach of Sec.2.1 we derive a bidirectional 
Poisson distribution for tunneling current, obtain a Schottky-like relation 
for the third moment and discuss its robustness. After that we discuss the 
relation of the counting statistics theory and the theory of photo-detection 
(Sec. 2. 3). Then we proceed to the problem of mesoscopic transport. In 
Sec. 3.1 we review the results on the DC transport and the derivation of 
a functional determinant formula for the counting distribution generating 
function. In Sec. 3. 3 we review the work on the AC transport statistics, and 
then consider the problem of mesoscopic pumping (Sec. 3. 4). The counting 
statistics for generic pumping strategy at weak pumping is given by a bidi- 
rectional Poisson distribution. We show that the Fano factor varies between 
and 1 as a function of the pumping fields phase difference. 

2. General approach 

The transmitted charge distribution can be characterized [18, 19] by elec- 
tron counting probabilities p n , usually accumulated in a generating function 1 



The function x(^) is 27r-periodic in the counting field A and has the property 
x(0) = 1 which follows from the probability normalization J2 n Pn = 1- The 
term "counting field" will be motivated in Sec.2.1, where a microscopic 
definition of x(A) is discussed in which A appears as a field that couples 
current fluctuations to a charge detector. 

The generating function (1) is particularly well suited for characterizing 
statistics of the distribution p n . The so-called irreducible correlators ((5n k )) 
(also known as cummulants) are expressed in terms of x(A) as 




(1) 




fc=i 




(2) 



The first two correlators in (2) give the mean and the variance: 



mi = n, 



771,2 = o~n 2 = n 2 — n 2 , 



(3) 



1 The function x(A) is also called a characteristic function [1]. 
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where f(n) stands for J2n f( n )Pn- The third correlator 2 

m 3 = ((5n 3 )} = h?= (n - nf (4) 

characterizes the asymmetry (or skewness) of the distribution p n . 

To illustrate the notion of a generating function, let us consider a Pois- 
son process. It describes charge transport at very low transmission, with 
uncorrelated transmission events. For the Poisson distribution 

p k = {f^l k[ m and X (A) = exp ( (e * - l)n) , (5) 

where n = It/ e is the average number of particles transmitted during time 
t, with / the time-averaged current and e the elementary charge. Comparing 
(5) with (2), one finds that all cummulants of the Poisson distribution are 
identical: ra k = n. 

Another useful example is binomial statistics. A binomial distribution 
arises when a fixed number N of independent attempts to transmit particles 
is made, each attempt successful or unsuccessful with probabilities p and 
q = 1 — p. The probability to transmit k particles in this case is determined 
by the combinatorial number C^- = N\/(N — k) \k\ of k successful outcomes. 
The probability distribution and the generating function in this case are 

Pk = C k N p k q N - k and X (A) = (pe lX + q) " . (6) 

The cummulants of the binomial distribution (6) can be found from (2) by 
expanding ln X in A: 

mi=pN, m2=pqN, m.3 = pq(q — p)N, ... (7) 

The binomial distribution (6) describes counting distribution of DC current 
noise for a single channel scatterer, such as point contact, at zero temper- 
ature (Sec. 3.1). 

Statistically independent processes result in a generating function given 
by a product of generating functions for constituting processes: x(A) = 
%i(A)...Xfe(A). For example, consider a biderectional Poisson distribution 
defined as a mixture of two independent Poisson processes transmitting 
particles in opposite directions with the rates n and n' . In this case, 

X2P(X) = exp ((e iA - l)n) • exp ((e~ iA - l)n') . (8) 

2 The relation between cummulants and correlators is generally more complicated than 
Eq.(4) for the third cummulant. For example, 7714 = ({Sq 4 )) = Sq 4 — 3 (5g 2 ) . 
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In Sec. 2. 2 we use the distribution (8) to describe statistics of tunneling 
current. In Sec. 3. 4 we show that it describes noise a mesoscopic pump. 

2.1. A MICROSCOPIC REPRESENTATION OF x(A) 

Here we discuss a microscopic definition of counting statistics for a physi- 
cal system. Adopting an inductive approach, we shall start with a specific 
model of current detector [20, 24]. We obtain the generating function %(A) 
for a particular current measurement scheme, and then argue that it de- 
scribes generic measurement. 

In a realistic noise measurement, e.g. in a mesoscopic wire or a point 
contact, the current fluctuations are not detected directly. Instead, the 
measurement is performed on the electromagnetic fluctuations (basically, 
voltage noise) induced by current fluctuations in the system. The electro- 
magnetic fluctuations have to be amplified before being detected. The con- 
version of underlying microscopic fluctuations due to fermions (electrons, 
fractional charges, etc.) into fluctuations of bosons (photons) is crucial, 
since Bose fields can be amplified without compromising noise statistics, 
while Fermi statistics is not consistent with amplification. 

Our goal is to clarify the microscopic picture of current fluctuations, 
rather than to describe realistic measurements. Thus we choose a gedanken 
measurement scheme well suited for that purpose. Consider a spin 1/2 
placed near an electron system and magnetically coupled to the electric 
current. We restrict the coupling to the spin z component, so that the 
system in the presence of the spin is described by Ti(q, p — ao^), where 
H(q,, p) is the electron Hamiltonian and a(r) is the spin vector potential 
scaled by e/c. 

The scheme of current detection using such spin dynamics can be moti- 
vated quasiclassically. A spin coupled to a time-dependent classical current 
I(t) by the interaction H = ^Xasl(t) will precess at the rate proportional 
to current, which turns the spin into an analog galvanometer. Indeed, if the 
spin-current coupling is turned on at t = 0, the spin will start precessing 
around the z axis with the precession angle 6(t) = A Jq I(t')dt' equal to the 
transmitted charge times A. The coupling constant A, so far arbitrary, will 
be associated with counting field below. 

In a fully quantum-mechanical problem, the spin evolution can be ob- 
tained from id = [a, H] . Since the spin-current coupling Hamiltonian com- 
mutes with (73, the spin dynamics can be found explicitly. For that we con- 
sider the transverse spin components a± = <j\ ± io<i and write the evolution 
equation i&± = [<j± , H] in the form 



i<j + = cr+P L W(q, p + a) - H(q, p - a)P T cr + , 
ia- = CT_P T W(q, p - a) - H(q, p + a)P[<J- , 



(9) 
(10) 
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with P^i = ^(1 + u 3) the up and down spin projectors. Here we used the 
raising/lowering properties of the operators a± and replaced 03 by 1 to the 
left of <t + and by —1 to the right of a + (and similarly for the <r_ equation). 

We consider a measurement which is performed during time interval 
< t' < t, i.e. start with a free spin, a. t <o = 0, couple it to the electron 
system at t = 0, maintain a finite coupling during < t' < t, and then turn 
it off. The expectation value of the transverse spin component at the time 
t, found by integrating Eqs. (9), (10), is 

(a+(t)) = (e iW{q ' p - a)t e ^ (q ' p+a) *)ei(^ + (0))s P in, (11) 

while (<7_(i)} = (a + (t)}*. Note that the result of the coupled spin and 
current evolution factors into a product of quantities that depend separately 
on electron dynamics and on the initial state of the spin, as indicated by 
the subscripts. 

The effect of current on spin precession is described by the dependence 
of the first term in Eq.(ll) on the gauge field a. To make contact with 
the quasiclassical discussion above, let us expand 7i in a, assuming it to 
be small. The result is H(q, p ± a) = W(q, p) ± aj, where j is electric 
current. Substituting this back in Eq.(ll) and passing to the interaction 
representation with respect to the Hamiltonian of fermions uncoupled from 
the spin, we rewrite the average (...} c i in (11) as 

^Texp (i £ J ai(t')d 3 r dt'^J Texp (-i J* J aj(t')d 3 r dt'^ . (12) 

Let us consider a specific form of a, taking it to be a pure gauge, § adl = 0, 
within the electron system, and nonzero near a particular surface (e.g. a 
5— function on the surface). For a classical current, ignoring noncommuta- 
tivity of current operators at different times, the expression (12) becomes 

(e~ im ) el with 0(t) = A J* I(t')dt' '. (13) 

Here I(t) is the total current through the surface and A = —2 J adl, where 
the integral is taken across the surface. The form of Eq.(13) agrees with 
what one expects for the precession phase factor averaged over classical 
current fluctuations. For n electrons transmitted through the system during 
the measurement time, the precession angle is 9(t) = An. This relates the 
average in (13) with the transmitted charge distribution: 

^ (t) ) el = £^> n . (14) 
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The relation with the spin precession in this case can be seen more clearly 
by combining the result (14) with Eq.(ll), 



and recalling the transformation rule a' + = e ie a + for spin rotation around 
the z-axis by an angle 0. This way of writing the result of spin evolution 
confirms the expected relationship between the charge counting probability 
distribution and the distribution of the spin precession angles. 

This discussion clarifies the meaning of the quantity (...) e i in Eq.(ll), 
linking it to the counting distribution generating function. Motivated by 
this, we use Eq.(ll) to give a microscopic definition of counting statistics. 
We rewrite the quantity (...) Keldysh partition function 



with the integral taken over the Keldysh time contour Co,t = [0 — > t — > 0], 
first forward and then backward in time. The counting field A is related 
to the gauge field a via A = Tf / adl, where the integration path goes 
across the region where scattering takes place and noise is generated (e.g. 
across the barrier in the point contact). The sign T indicates that the field 
a is antisymmetric on the upper and lower parts of the Keldysh contour. 
Because of that, even though a resembles in many ways an ordinary elec- 
tromagnetic gauge field (allowing for gauge transformations, etc.), it has no 
such meaning. We emphasize that a is really an auxiliary field describing 
coupling with a virtual measurement device, such as the spin 1/2 above. 

The microscopic formula (16), originating from the analysis of a coupling 
with spin 1/2, is in fact adequate for any ideal "passive charge detector" 
without internal dynamics. We shall use this formula below to obtain count- 
ing statistics for several physical situations of interest, including tunneling 
and mesoscopic transport. 

This still leaves some questions about universality and limitations of 
Eq.(16). Nazarov and Kindermann [37] considered a more general scheme 
of charge detection and recovered the expression (16). Although this is 
reassuring, Ref. [37] concludes that the detector back action is inevitable. 
Thus it is still desirable to study more realistic models of noise detection 
that include conversion of microscopic current noise into electromagnetic 
field (photons) as well as an amplifier. Some aspects of electron-to-photon 
noise conversion were studied by Beenakker and Schomerus [32]. 

It is also of interest to compare the back action effects in different mod- 
els. We argue that the above scheme is likely to describe noise measurement 




(15) 



n 




(16) 
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with the least back action, since coupling to the precessing spin 1/2 affects 
only the phase of electron forward scattering amplitude, without changing 
scattering probabilities. 

2.2. STATISTICS OF THE TUNNELING CURRENT 

The problem of the tunneling current noise provides a simple test for the 
microscopic formula (16). The starting point of our analysis will be the 
tunneling Hamiltonian H = H1 + H2 + V, where Hi^ describe the leads and 
V = J12 + J21 is the tunneling operator. The specific form of the operators 
Ji2-> J21 that describe tunneling of a quasiparticle between the leads will 
not be important for the most of our discussion. Both the discussion and 
the results for the tunneling current statistics obtained in this section are 
valid for a generic interacting many-body system. 

The counting field a in this case can be taken localized on the bar- 
rier, entering the Hamiltonian through the phase factors exp(±i / adX) = 
exp(±iA/2) of the operators J12, ^21- The tunneling operator then is 

V x = e^J 12 (t) + e-^)j 21 (t). (17) 

Here X(t) = ±A is antisymmetric on the Keldysh contour Cq^- 

In what follows we compute x(^) an d establish a relation with the Kubo 
theorem for the tunneling current [39]. For that, we perform the usual gauge 
transformation turning the bias voltage into the tunneling operator phase 
factor as J12 — > Ji2e~ teVt , J21 — ► J2ie teVt - Passing to the Keldysh interac- 
tion representation, we write 




Diagrammatically, the partition function (18) is a sum of linked cluster 
diagrams with appropriate combinatorial factors. To the lowest order in 
the tunneling operators J12, J21 we only need to consider linked clusters of 
the second order. This gives 

X (\)=e w ^\ ^^l/^j^WV)^'))^ ( 19 ) 

This result is correct for the measurement time t much larger than the 
correlation time in the contacts that determines the characteristic time 
separation t' — t" at which the correlator in (19) decays. 

There are several different contributions to the integral in (19), arising 
from t' and t" taken on the forward or backward parts of the contour Cqj ■ 
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Evaluating them separately, we obtain 

W(X) = f f {V-\V)Vtf'j) dt"dt> (20) 
- ff (Vx(t')Vx(t")) df'dt 1 - f ' ['(y-xit'W-xit'')) dfdt'. 

"0^0 " J t 

We substitute the expression (17) in Eq.(20) and average by pairing J\i 
with J21. This gives 

W(X) = (e iX -l)N^ 2 (t) + {e- iX -l)N 2 ^(t) (21) 

with 

N i^k = f f ' (J kj (t')Jjk(t")) dt'dt" = g jk t (22) 
Jo Jo 

the mean charge transmitted from the contact j to the contact k in a 
time t. Exponentiating Eq.(21) gives nothing but the bidirectional Poisson 
distribution X2p(A) defined by Eq.(8) with the transition rates given by 
n = Ni^ 2 = gut, n' = N 2 ^i = g 2 it, respectively. 

Eq.(8) yields interesting relations between different statistics of the dis- 
tribution. The cummulants ((Sn k )), obtained by expanding ln%2p(A) in A, 
are 

XX " I (512 +521 )t, keven. v ' 

Setting k = 1, 2 we express 312+521 through the time-averaged current and 
the low frequency noise spectral density 3 : 

312 - 921 = I/qo, 912 + 921 = S /2ql, (24) 

with qo the tunneling charge. Of special interest for us will be the third 
correlator of the transmitted charge 



((6q 3 )) = (6q-6q) . (25) 

For this correlator Eq. (23) gives ((Sq 3 )) = C%t with the coefficient C3 (the 
spectral density of the third correlator at uj = 0) related to current as 

C3 = ((5q 3 ))/t = qp. (26) 

3 The spectral density of the noise is defined through the symmetrized current correla- 
tor &s Sui — J { SI (t), 5/(0)} e wt dt. At w = 0, one can write So in terms of the variance 

of charge q(r) = f* +T I(t')dt' transmitted during a long time r as 5 = 7 <5? 2 (t). 
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We note that the relation (26) holds for the distribution (8) at any ratio 
(512 — 921)/ (912 + 521) of the mean transmitted charge to the variance. 

The quantities (22) have several general properties. First, by writing 
the expectation values (22) in a basis of exact microscopic states and using 
the detailed balance relation, we obtain 

N^z/N^i = 912/921 = exp( e y/fc B T) , (27) 

where V is the voltage applied to the contacts. Using this result to cal- 
culate the ratio of the first and second cummulants, Eq.(24), we have 
(512 — <72i)/(<7i2 + 521) = tanh(ey//cB^). This gives the noise-current re- 
lation 

S = 2q coth(eV/k B T) I (28) 

that holds for arbitrary eV/k^T. This relation was pointed out by Sukho- 
rukov and Loss [38] . 

Also, one can establish a relation of the quantities (22) with the Kubo 
theorem. We consider the tunneling current operator 

l(t) = -iq (Ji 2 (t) - J 2 i(t)) • (29) 

From the Kubo theorem for the tunneling current [39] , the mean integrated 
current ^{X{l/))d£ is 

qo f f ( [ J21 (0 J12 (<")] ) dt'dt" = q (N 1 ^ 2 -N 2 ^ 1 ) . (30) 

By writing Nj^k = 9jkt, we confirm the first relation (24). To obtain the 
second relation (24) we consider the variance of the charge transmitted in 
time t, given by ((5q 2 )} = ql ^{{i{t'),i{t")} + ) dt'dt" . This integral can 
be rewritten as 

f f {{j™W)> dt'dt" = N^+N^i , (31) 

which immediately leads to the second relation (24). 

We conclude that the tunneling current statistics, described by Eq.(8), 
are simpler than in a generic system. The current-noise relation, typically 
known in a generic system only at equilibrium (Nyquist) and in the fully 
out-of-equilibrium (Schottky) regimes, for the tunneling current is given by 
Eq.(28) at arbitrary eV/k B T. 

In contrast, the relation (26) obeyed by the third correlator (4) is com- 
pletely insensitive to the crossover between the Nyquist and Schottky noise 
regimes. The meaning of Eq.(26) is similar to that of the Schottky formula 
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So = 2((5q 2 }) = 2qoI. However, the Schottky current-noise relation is valid 
only when charge flow is unidirectional, i.e. at low temperatures k^T <C eV, 
since g\2j()2\ = exp(eV/fcBT), while Eq.(26) holds at any eV/k^T. 

In experiment, when the current-noise relation is used to determine 
the tunneling quasiparticle charge qo from the tunneling current noise, it 
is crucial to maintain low temperature k^T <C eV. The requirement of a 
cold sample at a relatively high bias voltage is the origin of a well known 
difficulty in the noise measurement. In contrast, the relation (26) is not 
constrained by any requirement on sample temperature. 

This property of the third moment, if confirmed experimentally, may 
prove to be quite useful for measuring quasiparticle charge. In particular, 
this applies to the situations when the I — V characteristic is strongly 
nonlinear, when it is usually difficult to unambiguously interpret the noise 
versus current dependence as a shot noise effect or as a result of thermal 
noise generated by non-linear conductance. This appears to be a completely 
general problem pertinent to any interacting system. Namely, in the systems 
such as Luttinger liquids, the I — V nonlinearities arise at eV > k-^T. 
However, it is exactly this voltage that has to be applied for measuring the 
shot noise in the Schottky regime. 

Finally, we note that the universality of the third moment is specific 
for the tunneling problem. In other situations, such as a point contact or a 
mesoscopic system, the third moment is temperature-dependent [34, 35, 36]. 

2.3. A RELATION TO THE THEORY OF PHOTODETECTION 

The statistics of tunneling particles was not specified in the above discus- 
sion, since everything said so far is good for both bosons and fermions. 
To illustrate this, here we discuss the relation of the present approach to 
the theory of photon counting [1, 2]. A system of photons interacting with 
atoms in a photon detector can be accounted for by a Hamiltonian of the 
form H = Hp + H a + V, where H p describes free electromagnetic field, H& 
is the Hamiltonian of atoms in the detector, and 



describes the interaction of photons with the atoms, i.e. the process of 
photon absorption and atom excitation. Here cik are the canonical Bose op- 
erators of photon modes, labeled by k, and bj are the operators describing 
excitation of the atoms. Since the operator V transfers excitations between 
the field and the atom systems, it can be interpreted as a "tunneling opera- 
tor." (The only difference is in the unidirectional character of the "current" 
induced by V, since photons can be only absorbed in the detector but not 




(32) 
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created.) This analogy allows one to use the formalism of Sec. 2.1 to study 
photon counting, and for that purpose we added a counting field in (32) 
(compare to Eq.(17)). 

Given all that, the generating function for photons has the form (18) 
which we rewrite to show an explicit dependence on the measurement time: 

Xt(A) = (uzl(t)U x (t)), U x {t) = Texp (-i j\{t')dt^j , (33) 

The task of evaluating the partition function (33) is simplified by the weak- 
ness of the photon-atom coupling. (Each atom is excited during the count- 
ing time t with a very small probability.) The expression (33) can thus be 
evaluated by taking into account the interaction of a photon with each of 
the atoms only to the lowest order. This is also similar to the tunneling 
problem. 

However, at this stage the similarity with tunneling ends, since photon 
coherence time can be much longer than the measurement time t. The 
method of Sec. 2. 2, based on the linked cluster expansion of lnx, should be 
modified to account for the long coherence times. Another complication is 
that the photon density matrix is not specified, since we are not limiting 
the discussion to thermal photon sources. 

We handle the partition function (33) by averaging over atoms, while 
keeping the photon variables free. As explained above, only pairwise aver- 
ages of atoms' operators are needed. We write them as (b^(t)bj'(t')) = 0, 

(bj(t)bj,(t')) = Tj5(t — t')5jjr, where Tj is a constant of the order of the 
excitation time of an atom, and the 5-function is actually a function of the 
width ~ Tj. (Typically, Tj is a very short, microscopic time.) 

Turning to the calculation, let us consider the difference x*+a(A) — Xt(A), 
with the time increment A large compared to Tj , but much smaller than the 
characteristic photon coherence time. Expanding U x {t + A) to the second 
order in A, we write it as 

l-ij t V x (t>)dt> - - jf jf V x {t')V x {t")dt'dt"\U x {t). (34) 

Substituting this in Eq.(33) along with a similar expression for Ull(t + A), 
and averaging over the atoms as described above, we obtain 

dtXt(X) = (Xt+A(A) - X t(A))/A = 5> fc (e a - 1) (uz l x (t)a\a k U x (t)) (35) 

k 

with rjk = J2j T j\ u j,k\ 2 t ne detector efficiency parameters. The solution of 
Eq.(35) has the form well known in optics [1, 2] : 

Xt(A) = II^ fe) ( A )' X? } (A) = ( : exp ( % i(e* A - l)a\a k ) : ) fc , (36) 

k 
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where : ... : is the normal ordering symbol and (...)& is the average over pho- 
ton density matrix. The product rule in (36) indicates that the counting 
distributions for different electromagnetic modes are statistically indepen- 
dent. The physical meaning of normal ordering is that each photon, after 
having been detected, is absorbed in the detector and destroyed. 

From Eq.(36), the counting probability of m photons in one mode is 



(k) 



Pm 



ml 



■(:(4a fc ) m e-" fc *°i afc :) fc . (37) 



Eqs.(36),(37) are central to the theory of photon counting [1]. Particularly 
interesting is the case of a coherent photon state |z), a|z) = z\z), with a 
complex z, corresponding to the radiation field of an ideal laser. In this 
case Eq.(37) yields the Poisson distribution p rn = e~ Jt (Jt) m /ml, J = i]\z\ 2 , 
which describes the so-called minimally bunched light sources. 



3. Counting statistics of mesoscopic transport 

Here we consider the problem of counting statistics in a mesoscopic trans- 
port. From now on we adopt the noninteracting particle approximation and 
use the scattering approach [40] , in which the system is characterized by a 
single particle scattering matrix. Depending on the nature of the problem, 
the matrix can be stationary or time-dependent. Even for noninteracting 
particles the problem of counting statistics remains nontrivial due to cor- 
relations between different particles arising from Fermi statistics. 

Counting statistics can be analyzed using the microscopic formula (16). 
However, there is a more efficient way of handling the noninteracting prob- 
lem. One can obtain a formula for the generating function x(^) i n terms 
of a functional determinant that involves the scattering matrix and the 
density matrix of reservoirs. Then for each particular problem one has to 
analyze and evaluate an appropriate determinant. Although functional de- 
terminants can be nontrivial to deal with, this approach is still much simpler 
than the one based directly on Eq.(16). 



3.1. STATISTICS OF THE DC TRANSPORT 

He we discuss the problem of time-independent scattering. We consider 
a conductor with m scattering channels describing states within one or 
several current leads. The scattering is elastic and will be characterized by 
a m x m matrix S. Although in applications so far the 2x2 matrices (i.e. 
the problems with two scattering channels) have been more common than 
larger matrices, the general determinant structure of x(A) will be revealed 
only for matrices of arbitrary size m. 



14 



L.S. LEVITOV 



For elastic scattering one can obtain x(A) from a quasiclassical argu- 
ment. In this case, particles with different energies contribute to counting 
statistics independently, and thus one can "symbolically" write 

x(A)=IIXe(A), i.e. X (A) = exp(t|ln Xe (A)^l), ( 38 ) 

where Xe(A) is the contribution of particles with energy e. The factor 2irh 
is written based on the quasiclassical phase space volume normalization, 
dV = dedt/2-irh. The quantity Xe(A) depends on the scattering matrix S 
and on the energy distribution raj(e) in the channels. 

To obtain x«(A) we introduce a vector of counting fields Xj, j = 1, m, 
one for each channel, and consider all possible multi-particle scattering 
processes at fixed energy. The processes can involve any number k < m of 
particles each coming out of one of the m channels and being scattered into 
another channel. Since the particles are indistinguishable fermions, no two 
particles can share an incoming or outgoing channel. One can then write 
Xe(A) as a sum over all different multiparticle scattering processes: 

Xe(A) = £ 6i^i+™+^-^--^)p il> ... >ifc | il> .. ., jk , (39) 
h,—,ik,ji,—,jk 

where the rate of k particles transition from channels i±, %k into channels 
ji, j k is given by 



Here S^''"'^ is an antisymmetrized product of k single particle amplitudes, 
which is nothing but the minor of the matrix S with rows ji,---,jk an d 
columns i±,...,ik- The product of rij and 1 — raj gives the probability to 
have k particles come out of the channels i±, ...,ik- 

An important insight in the structure of the expression (39) can be 
obtained by noting that it has a form of a determinant: 

Xe(\) = det (i - n £ + n e SZ\S x ) , (Sx)ij = e^'^S^. (41) 

Here h e is a diagonal m x m matrix of channel occupancy at the energy 
e and the counting field Xj enters in the phase factors of the matrix S\. 
To demonstrate that the expressions (39) and (41) are identical one has to 
expand the determinant (41) and go through a bit of matrix algebra. The 
formula (41) is particularly useful because, as we shall see below, it can be 
generalized to time-dependent problems. 
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Let us now focus on the simplest case m = 2 which describes transport in 
a point contact, with the two channels corresponding to current leads. The 
2x2 scattering matrix S contains reflection and transmission amplitudes. 
In this case, since there are only six terms in Eq.(39), the determinant 
formula (41) is not necessary. From Eq.(39) we obtain 



Xe(X) = (1 - m)(l - n 2 ) + (\S xl \ 2 + e^-^\S 2l \ 2 ) ni {l - n 2 ) 

+ {\S 22 \ 2 + e^ {Xl - X2) \S 12 \ 2 )n 2 {l - m) + |det S\ 2 nm 2 , (42) 



where the energy dependence of rij(e) is suppressed. By using the unitarity 
relations |Sn| 2 + |5 2i | 2 = 1, |detS| = 1, Eq.(42) can be simplified: 



Xe(X) = 1 +p(e tX - l)m(l - n 2 ) + p{e~ lX - l)n 2 (l - m) . (43) 



Here p = \S 2 \\ 2 = \Si 2 \ 2 is the transmission coefficient and A = X 2 — Ai. 
(We denote the transmission probability by p instead of a more traditional 
t to avoid confusion with the measurement time.) 

To obtain the full counting statistics integrated over all energies, one 
has to specify the energy distribution in the leads and use Eq.(38). We 
consider a barrier with energy-independent transmission and the leads at 
temperature T biased by voltage V. Then ni ;2 = np(e^f eV/2) with rap the 
Fermi function. 

At T = 0, since np(e) takes values and 1, for V > we have 



After doing the integral in Eq.(38) we obtain the binomial distribution (6), 
Xat(A) = (e iX p + 1 -p) N ®, with the number of attempts N(t) = eVt/2nh. 

This means that, in agreement with intuition, in the energy window eV 
the transport is just the single particle transmission and reflection, while 
the states with energies in the Fermi sea, populated in both reservoirs, are 
noiseless. (At V < the result is similar, with e tX replaced by e~ tX , which 
corresponds to the DC current sign reversal.) 

We note that the noninteger number of attempts N(t) = eVt/2 r K?i is 
an artifact of a quasiclassical calculation. In a more careful analysis the 
number of attempts is characterized by a narrow distribution Pjv peaked 
at N = N(t), and the generating function is a weighted sum ^2 N PnXn{X). 
Since the peak width is a sublinear function of the measurement time t 
(in fact, 5N 2 oc hit), the statistics to the leading order in t are correctly 
described by the binomial distribution. 

One can also consider the problem at arbitrary k^T jeV [20]. The inte- 
gral in Eq.(38), although less trivial, can be carried out, giving 




(44) 



X(A) = exp (-u + u_ N T ) , N T = t k B T/27rh, (45) 
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where 

u± = v ± cosh _1 (pcosh(« + iX) + (1 — p) coshu), v = eV/2ksT. (46) 

At low temperature k-gT <C eV, Eq.(45) reproduces the binomial statistics. 
At low voltage eV <C k^T (or high temperature) Eq.(45) gives the counting 
statistics of the equilibrium Nyquist noise: 

X (A) = e~ x * NT , sin(A*/2) = p 1 / 2 sin(A/2) . (47) 

Remarkably, even at equilibrium the noise is non-gaussian, except for spe- 
cial case of full transmission, p = 1, A* = A, when it becomes gaussian. 

3.2. STATISTICS OF TIME-DEPENDENT SCATTERING 

The time-dependent scattering problem describes photon-assisted trans- 
port. There are two groups of practically interesting problems: the AC- 
driven systems with static scattering potential, such as tunneling barriers 
or point contacts in the presence of a microwave field [17, 41], and the elec- 
tron pumps with time-dependent scattering potential controlled externally, 
e.g. by gate voltages [42]. 

Typically, the time of individual particle transit through the scattering 
region is much shorter than the period at which the system is driven. This 
situation is described, in the instantaneous scattering approximation, by 
a time-dependent scattering matrix S(t) that characterizes single particle 
scattering at time t. The question of interest is how Fermi statistics of 
many-body scattering states affects the counting statistics. 

One can construct a theory of counting statistics of time-dependent 
scattering [19] by generalizing the results of Sec. 3.1 for the statistics of a 
generic time-independent scattering. In particular, as we discuss below, 
the determinant formula (41), along with Eq.(38), allows a straightforward 
extension to time-dependent problems. In that, the generating function 
x(A) acquires a form of a functional determinant. 

Let us consider a scattering matrix S(t) varying periodically in time 
with frequency fi. The analysis is most simple in the frequency represen- 
tation [19], in which the scattering operator S has off-diagonal matrix ele- 
ments S u i^ with discrete frequency change oj' — to = nVt. In this approach 
the energy axis is divided into intervals nQ. < uo < (n + l)Q and each such 
interval is treated as a separate conduction channel. In doing so it is conve- 
nient (and in some cases necessary) to assign a separate counting field A n to 
each frequency channel, giving the counting field a frequency channel index 
in addition to the conduction channel dependence displayed in Eq.(41). 

Since the scattering operator conserves energy modulo multiple of h£l, 
the scattering can be viewed as elastic in the extended channel represen- 
tation, which allows to employ the method of Sec. 3.1. We note also that 
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the form of the determinant in Eq.(41) is not particularly sensitive to the 
size of the scattering matrix. Thus one can use it even when the number 
of channels is infinite, provided that the determinant remains well defined. 
This procedure brings (41) to the form of a determinant of an infinite size 
matrix. Determinant regularization can be accomplished by truncating this 
matrix at very high and low frequencies, thereby eliminating empty states 
and the states deep in the Fermi sea which do not contribute to noise. 

Finally, we note that the product rule (38) for x(A) is consistent with 
the determinant structure, since scattering processes at the energies dif- 
ferent modulo raO are decoupled. This allows to keep the answer for x(A) 
in the form of the determinant (41), where now the scattering operator S 
is considered in the entire frequency domain, rather than at the discrete 
frequencies w+nH. The resulting functional determinant has a simple form 
in the time representation: 

X (A)= det(l + n(t,t')(f x (t)-l)), f x (t) = sl x (t)S x (t), (48) 

where {S\)jji = e2^' _ V) Sjf as above, and h is the density matrix of 
reservoirs. The operator n, diagonal in the channel index, is given by 

n jf (t, t') = S jf J nj (huj) e^*'-') dw/2TT . (49) 

In general n(i, t') depends on the energy distribution parameters, such 
as temperature and chemical potential. In equilibrium, by taking Fourier 
transform of the Fermi function np(e — fi), one obtains 

e -ifj,j(t-t')/h 

ni(t,t') = — tt-, ; r— -, (3 = h/k B T. (50) 

n ' 1 2/3sinh(7r(t-t / + i < 5)//3)' ' ' y ' 

Finite bias voltage V s described by [i\ — ^2 = eV. At T = and V = 0, 
Eq.(50) gives n(t,t') = i/(2n(t - t' + id)). The result (48) holds for an 
arbitrary (even nonequilibrium) energy distribution in reservoirs. 

The functional determinant of an infinite matrix (48) should be handled 
carefully. One can show that, in a mathematical sense, the quantity (48) is 
well defined. For the states with energies deep in the Fermi sea, n = 1 and, 
since det (f A (t)) = 1 due to unitarity of S, these states do not contribute 
to the determinant (48). Similarly, since h = for the states with very 
high energy, these states also do not affect the determinant. Effectively, 
the determinant is controlled by a group of states near the Fermi level, in 
agreement with intuition about transport in a driven system. The absence 
of ultraviolet divergences allows one to go freely between different repre- 
sentations, e.g. to switch from the frequency domain to the time domain, 
which facilitates calculations [24, 21]. 
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The above derivation of the formula (48) based on a generalization of 
the result (41) for time-independent scattering might seem not entirely 
rigorous. A more mathematically sound derivation that starts directly from 
the microscopic expression (16) was proposed recently by Klich [43]. 

3.3. CASE STUDIES 

Here we briefly review the time-dependent scattering problems for which 
the counting statistics have been studied. From several examples for which 
x(A) has been obtained it appears that the problem does not allow a general 
solution. Instead, the problem can be handled only for suitably chosen form 
of the time dependence S(t). 

In Ref. [19] a two channel problem was considered with S(t) of the form 

Q< \= ( r t ' \ - ( B + be ~ iUT A _+ ^ inT 

b[T) -{t r' ) ~\A + ae- iQT -B-be mr 

which is unitary for any t provided |yl| 2 + \a\ 2 + \B\ 2 + \b\ 2 = 1, Aa + Bb = 0. 
The problem was solved by using the extended channel representation in 
the frequency domain, in which each frequency interval nVL < uj < (n + 1)S7 
is treated as a separate scattering channel, as discussed above. 

For the reservoirs at zero temperature and without bias voltage the 
charge distribution for m pumping cycles is described by 

X (A) = (l +Pi(e lA - 1) + P2 (e- a - I))™ (52) 

with pi = \a\ 4 /(\a\ 2 + |6| 2 ) and p 2 = |6| 4 /(|a| 2 + \b\ 2 ). This result means 
that at each pumping cycle one electron is pumped in one direction with 
probability p±, or in the opposite direction with probability p 2 , or no charge 
is pumped with probability 1 — p\ — p 2 . The multiplicative dependence of 
x(A) on the number of pumping cycles m indicates that the outcomes of 
different cycles are statistically independent. One can thus view (52) as a 
generalization of the binomial distribution (6). 

The problem (51) was also studied in Ref. [19] at a finite bias voltage, 
when the counting distribution is not as simple as (52). To describe the 
result, for a given bias voltage V we find an integer n such that nf < 
fV < (n + 1)/, where / = Q/2ir is the cyclic frequency in (51). Then for 
a long measurement time t S> O -1 the counting distribution is 

X(A)=X^ > (A)-X^ 1 (A), (53) 

where 7V< = (%V-nf) t, iV> = ((n + l)/-f V) t, and the functions Xn(A) 
are finite degree polynomials in e ±lX . The form of Xn(A) depends on A, B, 
a, and b (we refer to Ref. [19] for details). 
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The product rule (53) means that the cummulants of the distribution 
x(A) depend on V in a piecewise linear way, m^V) = N > m^ + N < m^ l+1 \ 
with cusp-like singularities at eV = nhf = nMl. These singularities are 
generic for the noise in photo-assisted phase-coherent transport [44, 17, 41]. 

Another time-dependent problem for which solution can be obtained in 
a closed form is mesoscopic transport in the presence of an AC voltage [24] . 
The scatterer in this case is a time-independent 2x2 matrix, while the 
voltage V{t) enters in the phase factors of the density matrix in (48): 

ni, 2 (t, = e±i W)-rt»n<$(t, t f ), <p(t) = %V(t) (54) 

(compare this with the formula (50) for constant bias voltage). 

The counting distribution (48) for a family of such problems has been 
studied in Ref. [21]. It was noted earlier [45] that noise is minimized at fixed 
transmitted charge for a special form of time-dependent voltage: 

Each of the Lorentzian voltage pulses (55) corresponds to a 2tt phase change 
in tp(t). Interestingly, the noise-minimizing pulses (55) have large degener- 
acy: they produce noise which is insensitive to the pulses' widths and 
peak positions t^. This calls for an interpretation of the pulses (55) as 
independent attempts to transmit charge. Not surprisingly, the counting 
statistics for such pulses was found to be binomial: 

X (A) = (l + t(e iA -l)) m (56) 

with t the transmission constant. The lowest possible noise for a current 
pumped by voltage pulses is thus equal to that of a DC current with the 
same transmitted charge. 

The method of Ref. [21] also allows to find the distribution for an ar- 
bitrary sum of the pulses (55) with alternating signs. For example, two 
opposite pulses 

h ( 2n 2r 2 



give rise to the counting distribution 

X(A) = 1 - 2F + F(e iX + e~ iX ) , F = t(l - t)\ Z J ~_ ^ | 2 (58) 

with zip = hp + irip. The quantity |...| 2 is a measure of pulses' overlap in 
time, varying between for a full overlap and 1 for no overlap. For nonover- 
lapping pulses, %(A) factors into (t e tX + 1 — t){te~ lX + 1 — t), in agreement 
with the interpretation of a binomial distribution for independent attempts. 
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3.4. MESOSCOPIC PUMPING 

A DC current in a mesoscopic system, such as an open quantum dot, can 
be induced by pumping, i.e. by modulating its area, shape, or other param- 
eters [46, 48, 47]. After pumping was demonstrated experimentally [42], it 
came into the focus of mesoscopic literature (for references see [47, 49]). In 
particular, Brouwer made an interesting observation that the time-averaged 
pumped current is a purely geometric property of the path in the scattering 
matrix parameter space, insensitive to path parameterization. 

Transport through a mesoscopic system is described [40] by a scatter- 
ing matrix S which depends on externally driven parameters and varies 
cyclically with time. The matrix S(t) defines a path in the space of all 
scattering matrices. For a system with m scattering channels, the matrix 
space is the group U(m) = SU(m) x U(l). In an experiment one can, in 
principle, realize any path in the space of scattering matrices. 

Counting distribution for a parametrically driven open system was dis- 
cussed by Andreev and Kamenev who adapted the results [19] 

obtained for specific pumping cycles [Eqs.(51),(52)]. However, since the 
relation between the path in the scattering matrix space and the external 
pumping parameters is generally unknown, only the results valid for generic 
paths are of interest in this problem. 

Here we consider the weak pumping regime, when the path S(t) is a 
sufficiently small, but otherwise arbitrary loop, and show that in this case 
the counting distribution is universal [22] , taking the form of bidirectional 
Poisson distribution (8). From that, we obtain the dependence of the noise 
on the amplitude and relative phase of the voltages driving the pump. 

Before turning to the calculation, we discuss general dependence of 
counting statistics on the path in matrix space. Different paths S(t), in 
principle, give rise to different current and noise. However, there is a re- 
markable property of invariance with respect to group shifts. Any two paths, 

S(t) and S'(t) = S(t)S , (59) 

where So is a time- independent matrix in U (m) , give rise to the same count- 
ing statistics at zero temperature. We note that only the right shifts of the 
form (59) leave counting statistics invariant, whereas the left shifts gen- 
erally change it. One can explain the result (59) qualitatively as follows. 
The change of scattering matrix, S(t) — ► S'(t) = S(t)So, is equivalent to 
replacing states in the incoming scattering channels by their superpositions 
tp a = SQpipP . At zero temperature, however, Fermi reservoirs are noiseless 
and also such are any their superpositions. Correlation between superpo- 
sition states of noiseless reservoirs is negligible, while current fluctuations 
arise only due to the time-dependent scattering. Therefore, noise statistics 
remain unchanged. A simple formal proof of the result (59) is given below. 
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For a weak pumping field it is sufficient to evaluate (48) in the time 
domain by expanding lndet(...) in powers of 5S and keeping non- vanishing 
terms of lowest order. In doing so, however, we preserve full functional 
dependence on A which gives all moments of counting statistics. We write 
S(t) = e A ^S^ with antihermitian A(t) representing small perturbation, 
tr^^4 <C 1. Here is scattering matrix of the system in the absence of 
pumping. Substituting this into (48) one obtains 

f x ( t ) ee if + 5T x (t) = S^e- A -^e A ^S x 0) (60) 

with fj 0) = S m x f S ( x 0) and A\(t) = e^ 3 A{t)e~ i ^ <n (here a 3 equals +1 for 
the left and —1 for the right channel). Now, we expand (48): 

lnx(A) = IndetQo + tiR - \trR 2 + ^-tri? 3 - ... , (61) 

where Q$ = 1 + h(T x ^ — 1) and R = QQ l fi8T x . At zero temperature, from 
n 2 = h it follows that detQ = 1 and R = S ( x 0) ~ l h (e' A -^e A ^ - l) S { x ] . 
Therefore, 

lnx(A) = tr hM - -tr(nM) 2 + -tr(nM) 3 - ... , (62) 

where M = e~ A - x ^e Ax ^ — l. Note that at this stage there is no dependence 
left on the constant matrix S^°\ which proves the invariance under the 
group shifts (59). 

We need to expand (62) in powers of the pumping field, which amounts 
to taking the lowest order terms of the expansion in powers of the matrix 
A{t). One can check that the two 0{A) terms arising from the first term in 
(62) vanish. The 0(A 2 ) terms arise from the first and second term in (62) 
and have the form 

ln X = ^tr (n (a 2 _ x + A 2 -2A_ X A X )) - hr(hB x ) 2 (63) 

with B\(t) = A\(t) — A_\(t). At zero temperature, by using h 2 = h, one 
can bring (63) to the form 

X - tr (n [Ax, A_ x ]) + \ (tr (n 2 B 2 ) - tr(hB x ) 2 ) . (64) 

The first term of (64) has to be regularized in the Schwinger anomaly 
fashion, by splitting points, t',t" = t ± r//2, which gives 

y n (t',t")tr (A_ x (f)A x (t f ) - A x {t")A_ x {t')) dt . (65) 
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Averaging over small r\ can be achieved either by inserting in (65) additional 
integrals over t', t", or simply by replacing A\(t) — > \ (A\(t) + A\(t')), etc. 
After taking the limit rj — > 0, Eq.(65) becomes 



i 

8tt. 



■ ^ tr (A-\dtA\ — A\dtA-\) dt . (66) 

The second term of (64) can be written as 

1 r rtr(B x (t)-B x (t')) 2 



4(2vr) 2 II (t - f ) 



-dtdt' . (67) 



We decompose A = clq + z+z^ with respect to the right and left channels, so 



that [03, ao] = 0, [a 3 ,z] = -2z, 



2zt. Then A x = e'^Ae^ 3 



a + e 1 2 z^ + e l z z, B x = (e 1 ? — e W, W = z^ — z. Substituting this 
into (66) and (67) one rewrites Eqs.(66),(67) in terms of W(t): A 



sin A 

and 



j>tv{[a^W]d t W)dt (68) 



(1-cosA) ft tr (W{t)-W{t'))\ u ^, 

j j — — d d ■ m 

From that we obtain the counting distribution for one pumping cycle: 

X (A) = exp (u(e lA - 1) + v{e~' x - 1)) (70) 

with the transmitted charge average / = e(u— v ) and variance J = e 2 (u+v ), 
related to the noise spectral density by <So = Jfl/ir (see Sec. 2. 2). 

The parameters u and v in (70) can be expressed through z(t) and z\t) 
as follows. Let us write z{t) as z + (t) + z-(t), where z + {t) and z_(i) contain 
only positive or negative Fourier harmonics, respectively. Then 



1 

u = — 

47T 



I 

4tt 



j> tr (zLd t z+ - z+d t zVj dt = ^ oj trz^z^ , (71) 
I tr (z-dtz+ ~ z\d t z^j dt = -J2^ ^_ w z w . (72) 



oj<0 



This demonstartes that u > and v > 0. It is straightforward to show that 
(68) equals isinA(n — v), whereas (69) equals (cos A — l)(u + v). 

4 Eq.(68) is essentially identical to the result obtained by Brouwer for the average 
pumped current [48]. The integral in (68) is invariant under reparameterization, and thus 
has a purely geometric character determined by the contour S(t) in U(m). 
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Now we consider a single channel pump, S(t) G U{2). In this case, 
z = z(t)a- and = z*(t)a+. For a harmonic driving signal, without loss 
of generality, one can write 

z (t) = Zl Vi cos(fit + 6) + z 2 V 2 cos(fit), (73) 

where Vi j2 are pumping signal amplitudes, and complex parameters zi j2 
depend on microscopic details. From (71) we find the rates 



1 

U= 4 



z x V x e Ld + z 2 V 2 



z x V x e- id + z 2 V 2 



(74) 



Both u and v can vanish at a particular amplitude ratio V±/V 2 and phase 
9. When this happens, the two Poisson processes (70) are reduced to one, 
and the current-to-noise ratio gives elementary charge, I/J = ie" 1 . This 
happens at the extrema of I/J as a function of w = {Vi/V 2 )e l6 , for (74) 
reached at w = —z 2 /zi, —z 2 /z±. 

Reducing the counting statistics (70) to purely poissonian by varying 
pumping parameters is possible, in principle, for any number of channels n. 
However, since the number of parameters to be tuned is 2n 2 , this method is 
practical perhaps only for small channel numbers. Although the method is 
demonstrated for non- interacting fermions, we argue that it can be applied 
to interacting systems as well. Poisson statistics results from the absence 
of correlations of transmitted particles, which must be the case in a generic 
system, interacting or noninteracting, at small pumping current. Using the 
dependence of the rates u, v on the driving signal to maximize I/J, i.e. to 
eliminate one of the two Poisson processes (70), one could then obtain the 
charge quantum in the standard way as e = J/ 1. 

To summarize the results of this section, in the weak pumping regime 
the distribution of charge transmitted per cycle is of bidirectional Poisson 
form (8), i.e., it is fully characterized by only two parameters, average 
current and noise. The current to noise ratio I/J, scaled by elementary 
charge, varies between 1 and —1, depending on the relation between driving 
signals phases and amplitudes. Thus the quantity max(|/|/J) gives the 
inverse of elementary charge without any fitting parameters. Polianski et 
al. [49] recently studied the dependence of I / J on the mesoscopic scattering 
ensemble parameters, and found that, within the random matrix theory, the 
nearly extremal values close to ±1 can be reached with finite probability. 
This may permit to use the noise in a pump to measure quasiparticle charge 
in open systems, such as Luttinger liquids. 
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